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A  TWO-PARAMETER  CLASS  OF  BESSEL  WEIGHTINGS  WITH  CONTROLLABLE 
SIDELOBE  BEHAVIOR  FOR  LINEAR,  PLANAR-CIRCULAR,  AND  VOLUMETRIC- 
SPHERICAL  ARRAYS;  THE  IDEAL  WEIGHTING-PATTERN  PAIRS 

INTRODUCTION 

A  wide  variety  of  time-domain  weightings  for  spectral  analysis,  whose 
frequency-domain  windows  have  very  good  sidelobe  behavior,  have  been  presented 
in  [I,2J.  Since  the  basic  mathematic s  descrioing  the  response  of  a  weignted 
linear  array  can  also  be  written  as  a  Fourier  transform,  these 
weighting -window  pairs  have  immediate  application  to  one-dimensional  array 
processing  as  well  as  spectral  analysis. 

Most  of  the  weighting-window  pairs  in  [1,2]  have  no  parameters  in  their 
design  equations;  that  is,  the  windows  are  fixed  and  cannot  be  altered,  as  for 
example,  in  the  Hanning  and  Hamming  windows.  A  few  windows,  such  as  the 
Dolph-Chebyshev  and  Kai ser-Bessel  [3,4],  do  have  a  single  parameter  in  their 
design  equations  that  allows  for  a  tradeoff  between  the  mainlobe  width  and  the 
ratio  of  mainlobe  to  peak  sidelobe.  However,  neither  have  any  control  over 
the  rate  of  decay  of  the  sidelobes,  the  Dolph-Chebyshev  case  having  no  decay, 
and  the  Kai  ser-Bessel  case  a  6  dB/octave  decay.  It  is  obvious  that  in  order 
to  control  both  the  sidelobe  decay  and  the  mainlobe-to-peak-sidelobe  ratio,  a 
two-parameter  family  of  weightings  is  necessary.  And  it  is  desirable 
(although  not  necessary)  for  the  window  to  possess  a  simple  analytical  form 
that  can  be  easily  understood  and  evaluated  for  a  range  of  parameter  values. 
Such  a  class  of  Bessel  weightings  is  presented  in  this  report. 

For  the  array  application,  the  weighting  is  applied  as  a  multiplicative 
factor  in  the  spatial  domain;  the  response  to  plane  wave  arrivals  from  various 
directions  is  called  the  pattern,  rather  than  the  window.  Here  we  will  give  a 
two-parameter  family  of  weighting-pattern  pairs  for  use  with  arrays  in  one, 
two,  or  three  dimensions,  and  shall  indicate  the  ideal  weightings  and 
corresponding  patterns  in  all  cases.  Special  cases  of  this  family  will  be 
shown  to  include  some  of  the  weightings  that  are  currently  employed  in  array 
and  signal  processing. 
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RESPONSES  OF  ONE-,  TWO-,  AND  THREE-DIMENSIONAL  ARRAYS 


LINEAR  ARRAY 

We  consider  a  continuous  line  array  located  on  the  x-axis  in  the  range 
( -R,R)  and  subject  to  symmetric  weighting  w1( x )  f0r  |xl  <  R.  For  a 
single-frequency  plane  wave  of  wavelength  x,  arriving  at  angle  da  relative 
to  the  normal  to  the  line  array,  the  array  voltage  response,  by  use  of 
time-delay  steering  to  look-angle  ,  is 

g{u)  =  J'  dx  exp  ||i2irj(  sinila-sin^  )j  w^x).  (1) 

-R 

By  letting  s  =  x/R  and  by  using  the  symmetry  of  weighting  wlt  we  Can  express 
response  (1)  as 


where  normalized  weighting 


w(  s)  -  (2w/*R  Wj(Rs), 
and  dimensionless  parameter 


(2) 


(3) 


u  =  2i3(  si  nd  -  si  ndu )  ( 4 ) 

X  3  a 

incorporates  the  relevant  features  of  array  geometry,  look  angle,  and  the 
arrival  wavelength  and  angle. 

Thus  the  response  (2)  of  a  line  array  is  a  cosine  transform  of  the 
normalized  weighting.  As  an  example,  rectangular  weighting  yields  response 
g(u)  prooortional  to  sin  u/u,  which  has  its  first  few  nulls  at  u  *  *,  2*,  3w. 
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PLANAR  ARRAY 

The  voltage  re^ionse  of  a  continuous  planar-circular  array  of  radius  R,  to 
a  plane  wave  of  wavelength  x  arriving  at  (polar,  azimuthal)  angles  (da>  ga) t 
and  subject  to  weighting  which  depends  only  on  the  distance  from  the  center  of 
the  array,  is  derived  in  appendix  A,  culminating  in  the  result  (A-ll).  It  is 


1 

g(u)  ■  X  ds  s  J  (us)  w(  s), 
0  0 

where  w ( s)  is  the  noma lized  weighting  and 


(5) 


u  5  2*5.  j^sin2^  +  sin2da-  2sin^sinda  cos^-ea^j  .  (6) 

Here  (dj,  )  are  the  (polar,  azimuthal)  look  angles;  that  is,  the  response 
(5)  of  a  planar-circular  array  is  a  zero-th  order  Bessel  transform  of  the 
normalized  weighting.  As  an  example,  rectangular  weighting  w  yields  response 
g(u)  proportional  to  Ji(u)/u,  which  has  its  first  few  nulls  at  u  5  3.83, 

7.02,  10.17. 


VOLUMETRIC  ARRAY 


The  voltage  response  of  a  continuous  volumetric- spherical  array  of  radius 
R,  with  weighting  dependent  only  on  the  distance  from  the  center  of  the  array, 
is  derived  in  appendix  B.  The  result  is  given  by  (B-10)  in  the  form 

9(u)  -  Ms),  (7, 


where  now 


h  2*5 


6 


2-2co 


s^cos^ 


-  2sin 


«fesinda 


(8) 


The  other  parameters  are  as  explained  in  the  previous  subsection. 
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Equation  (7)  has  the  basic  form  of  a  sine  transform.  As  an  example, 
rectangular  weighting  w  yields  response  g(u)  proportional  to 

sin  u-u  cos  u  ln 

3  , 

u 

which  has  its  first  few  nulls  at  u  =  4.49,  7.73,  10.90. 

UNIFIED  FORM  FOR  ARRAY  RESPONSES  IN  DIFFERENT  DIMENSIONS 

The  results  in  (2),  (5),  and  (7)  for  the  array  voltage  response  in  one, 
two,  and  three  dimensions,  .reflectively,  appear  to  be  quite  different. 
However,  they  can  all  be  written  in  the  basic  form 


* 

g(u)  =  |  ds  K(u,s)  w(s), 


where  the  kernel 


K(u,s)  =  sgj  J  (us)  = 


U  1 

H  cos(us)  for  a  »  -  j 
s  J0(us)  for  u  =0 

fr.  ■  + 


Here  Jw  is  a  Bessel  function  of  order  u,  and  we  have  used  [5;  10.1.1, 

10.1.11,  10.1.12],  Thus  all  the  responses  are  basically  Bessel  transforms  of 
the  normalized  weighting  w,  with  the  correspondences  given  in  the  following 
table. 

Table  1.  Identification  of  Values  of  y  in  (11) 


Number  of  Dimensions  Value  of  u 


4 


1 

2 

3 


-1/2 

0 

1/2 


TR  6761 


If  we  substitute  (11)  in  (10),  we  have  the  explicit  result  for  the 
re  sponse  pattern: 


1 

g(u)  a  J1  ds  s^jj  J(us)  w(s). 
0 


(12) 


Inspection  of  the  properties  of  the  Bessel  function  reveals  that  g(u)  as  given 
by  (12)  is  even  in  u;  see  [5;  9.1.10].  Thus  we  only  need  to  investigate  g(u) 
for  u  >  0. 
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HAMKEL  TRANSFORM  PAIRS 

The  Hankel  transform  pair  of  order  u  is  given  by  [6;  page  136] 

00 

f(u)  =  (  ds  (us)^J  (us)  F ( s)  for  u  >  0, 

v  U 

0 
<30 

F(s)  =  j*  ds  (  su)*  J(  su)  f(u)  for  s  >  0.  (13) 

0 

Thus  knowledqe  of  either  function  f  or  F  for  positive  arguments  enables 
determination  of  the  other  function  by  an  integral  transform.  Under  the 
identification 

F(s)  =  s  Ms),  f(u)  *  u  g (u),  (14) 


(13)  takes  the  form 


g(u) 

w(s) 


CO 


us)  M  s)  for  u  > 


0, 


( su)  g{ u)  for  s  >  0. 


(15) 

(16) 


Equation  (15)  is  more  general  than  (12)  in  that  it  allows  for  weighting 
Ms)  to  be  nonzero  for  s  >  1.  Equation  (16)  is  complementary  in  that,  given  a 
desired  pattern  g(u),  it  indicates  what  weighting  Ms)  is  required  for  s  >  0. 
However,  if  we  attempt  to  specify  some  desirable  pattern  g(u),  and  then  solve 
(16)  for  the  required  weighting  M  s) ,  it  will  generally  turn  out  that  the 
resultant  Ms)  will  be  nonzero  for  s  >  1.  Thus  not  any  pattern  g(u)  can  be 
selected  if  we  Insist  on  a  finite- support  weighting  Ms);  rather,  desirable 
candidate  patterns  can  be  substituted  In  (16)  and  the  corresponding  weighting 
Ms)  evaluated  to  see  if  it  is  zero  for  s  >  1.  If  not,  the  candidate  pattern 
is  disallowed  and  must  be  modified  or  discarded?  We  will  use  precisely  this 


procedure  in  a  later  section  when  we  determine  the  weightings  that  realize  the 


ideal  patterns  in  various  numbers  of  dimensions. 


if  pattern  g(u)  yields  w(  s) 
a  modified  weighting  a411  CM 


*  0  for  s  >  a,  the  scaled  pattern  g(u/a)  yields 
as),  which  is  zero  for  s  >  1,  as  desired. 
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DEFINITION  OF  TWO  BESSEL  FUNCTION  RATIOS 

It  will  be  very  convenient  notationally  in  the  following  to  employ  the 
shorthand  notations 


JU) 


co  2  ^ 

_1  ^  (-z  /4) 

2“  kiHa+l+k) 


( z  /4) 
k'.TU+l+k) 


(17) 

(18) 


for  these  two  Bessel  function  ratios;  these  types  of  functions  have  already 
been  encountered  in  (11),  (12),  (15),  (16).  They  are  extensions  of  the 
^(z)  functions  di  scussed  in  [7;  page  56].  Both  functions,  ^(z)  and 
ddcU),  are  single-valued  and  are  entire  in  z  for  any  a,  as  well  as  being 
entire  in  a  for  any  z  [5;  6.1.3,  9.1.1].  Special  cases  of  these  functions 
that  are  useful  here  are  given  by  [5;  9.6.6,  10.1.1,  10.1.11,  10.1.12, 
10.2.13,  10.2.14] 


l#a(0)  =  Jk(0)  =  — — -  , 

“  Va  2aP(  a+1 ) 

(J0(z)  =  Iq(z),  dd^z)  =  —z—  ,t^_1( z)  =  zl^z), 

Ki <«>  -(if—  j-  s1nhz  •  i-vt*  -gf<«  * 


(3+z2)  sinh  z  -  3z  cosh  z 


nh  z  -  cosh  z) 


2(z)  *(-/  [(3+z  )  cosh  z  -  3z  sinh  z]. 


A  useful  prooerty,  vrfiich  is  obvious  from  (17)  and  (18)  and  which  will  find 
frequent  application  here,  is 
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J)  (±1z)  -(J?(Z),  d  (±12)  =5-  (z).  (20) 

Go.  a  a  (/a 

These  relations  hold  true  for  all  values  of  z,  real  or  complex.  Properties 
similar  to  (19)  can  be  obtained  for  the  j^a(z)  functions;  for  example, 

J--l/2<2>  -(ij cosMtfil  =(v)  cos(z>. 

Other  useful  properties  follow  from  the  use  of  [5;  9.1.30  and  9.6.28], 


j(<z)  •  Kiz)‘z^iu)- 

and  from  [5;  9.1.27  and  9.6.26], 

iljz)  "  -ifjUu,  -  2(a-l)(J!a_1U>J 
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A  CLASS  OF  8ESSEL  WEIGHTINGS 


The  class  of  3essel  weightings  that  we  are  suggesting,  regardless  of  the 
number  of  dimensions  of  the  array.  Is  given  by 


w(  s) 


for  0  <  s  <  1,  where  v  >  -1,  B  >  0, 


(22) 


and  w(  s)  =  0  for  s  >  1.  There  are  two  real  parameters  in  this  class,  namely, 
v  and  B.  For  the  fiecial  case  of  v  =  0,  this  weighting  is  already  known  as 
Kai ser-Be ssel  [3,4];  its  pattern  has  nearly  the  optimun  energy  content  within 
a  specified  bandwidth  in  the  one-dimensional  application  to  spectral  analysis. 


Substitution  of  (22)  in  (12)  or  (15)  yields  the  closed  form  for  the 
pattern  [6;  page  30,  fourth  integral,  and  5;  9.6.3] 


Here  we  have  employed  definitions  (17)  and  (18)  and  used  (20).  The  condition 
on  u  guarantees  convergence  of  integral  (12)  at  s  =  0,  vWiereas  the  condition 
on  v  guarantees  convergence  of  integral  (12)  at  s  *  1.  The  first  form  of 
(23B)  is  more  convenient  computationally  for  u  2  B,  whereas  the  second  form  i s 
more  convenient  for  B  >  u. 


Weighting-pattern  pair  (22)-(23)  are  the  fundamental  results  for  the  class 
of  Bessel  weightings  under  consideration.  They  apply  to  one-,  two-,  or 


three-dimensional  arrays  when  u  Is  specialized  to  -  -j,  0,  or  ■£,  reflectively, 
and  when  u  is  interpreted  as  (4),  (6),  or  (8),  reflectively.  The  parameter  B 
i s  nonnegative  and  will  be  shown  to  control  the  ratio  of  mainlobe  to  peak 
si  delobe.  For  the  case  of  one  dimension,  u  *-£,  the  kernel  of  transform  (23A) 
1  s  a  cosine  (  see  the  top  line  of  (11) );  in  thi  s  special  case, 


the  pattern^^ 
by  Roy  Streit  [8]. 


was  also  Independently  and  simultaneously  discovered 


TR  6761 


WEIGHTING  CHARACTERISTICS 

Weighting  (22)  is  positive  for  0  S  s  <  1  since  v  >  -1  and  B  £  0;  see 
[5;  9.6.10].  In  addition,  it  is  zero  and  therefore  continuous  at  s  =  1  if 
v  >  0.  In  fact  [5;  9.6.7], 


w(  s)~ 


(l-s)v 

Hv+1) 


as  s  — »  1- . 


(24) 


For  v  >  0,  weighting  w(  s)  i  s  monotonic  ally  decreasing  in  s  on  (0,1);  see 
appendix  C.  However,  for  -1  <  \>  <  0,  w(  s)  possesses  an  integrable  singularity 
at  s  3  1. 


Examples  of  weighting  (22)  are  plotted  in  figures  1-6  for  v  =  2,  1.5,  1, 
.5,  0,  -.5,  respectively.  Figure  5,  for  v  3  0,  corresponds  to  the 
Kai ser-Be ssel  weighting  [3].  For  the  larger  values  of  v,  the  weighting  blends 
smoothly  to  zero  at  s  =  1,  but  for  the  smaller  values  of  v,  the  behavior  of 
w(  s)  i  s  more  irregular  at  s  3  1,  being  discontinuous  for  v  =  0  and  infinite 
for  v  <  0.  The  larger  values  of  B  lead  to  smoother  functions  that  are 
Gaussian-like;  in  fact,  for  s  <  1  [5;  9.7.1], 


A  s)~ 


as  B-* * . 


More  generally. 


w(  s)~  (2w)4B_2v(bVi-s2'/  exp^fl-s2)  as  B^l- s2’-* «. 
The  opposite  limit  for  small  B  is 


w(  s) 


for  B  3  0. 


(25) 


(26) 


(27) 


10 


s 

Figure  5.  Normalized  Weighting  for 


V  =  -.5 


s 


.6 


:ed  Weighting  for  v 
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RESPONSE  fATTERN  CHARACTERISTICS 
The  pattern  was  given  by  (23)  as 

<»> 

where  the  composite  order  of  the  Bessel  function  1  s 

a=ii  +  v  +  l.  (29) 

The  asymptotic  behavior  of  the  pattern  (28)  for  large  u  is  available  from  (17) 
and  [5;  9.2.1]: 


g(u)~ 


as  u-*«°. 


(30) 


Since  g  is  proportional  to  the  array  voltage  response,  (30)  corresponds  to  a 

decay -v  6a  +  3  dB/octave  as  u  -*•  <»•  (31) 


Expressed  in  terms  of  the  original  dimension-parameter  n  and  weighting 
parameter  v,  this  i  s,  from  (29), 


decay  —  5u  +  6v  +  9  dB/octave  as  u  — *  00 
[6v  +  6  dB/octave  for  one  dimension 
‘  \6v  +  9  dB/octave  for  two  dimensions 
(6v  +  12  dB/octave  for  three  dimensionSj 


(32) 


Thus  the  greater  the  number  of  dimensions,  the  faster  is  the  rate  of  decay  of 
the  sidelobes  of  the  response,  for  a  common  value  of  weighting-parameter  v. 
Each  additional  dimension  adds  a  3  dB/octave  decay,  for  a  fixed  v. 
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Special  cases  of  the  pattern  (28)  are  available  upon  use  of  (19)  and  (20); 
for  example. 


g(u)  = 


for  a  =  0 


for  <x  =  ^ 


for  a  =  1 


(33) 


All  of  these  relations  are  valid  for  all  u,  whether  u  is  larger  or  smaller 
than  B;  of  course,  the  former  form  in  each  line  i s  more  convenient  for  u  <  B, 
while  the  latter  1 s  more  convenient  for  u  2  3.  The  third  result  in  (33),  for 


a  =  -j,  includes  the  pattern  in  one  dimension  (u  s  v  =0)  for  the 
Kai  ser-Bessel  weighting  10(q^i-s2)  for  0  <  s  <  1. 


The  special  case  of  a  s  -  j  in  (33)  deserves  extra  attention;  this  case 
will  be  called  the  ideal  pattern: 


(34) 


The  plot  in  figure  7  reveals  that  the  sidelobes  are  all  equal,  and  that  the 
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Figure  7.  Ideal  Pattern  gi(u) 


voltage  ratio 

mainlobe  level  _ 

side  lobe  level  ~  C0Sh(B)’  (35) 

The  mainlobe  width,  as  measured  to  the  point  vrfiere  the  mainlobe  first  decays 
to  the  eventual  sidelobe  level,  is 

mainlobe  width  =  B.  (36) 

The  abscissa  u  Is  given  by  (4),  (6),  or  (8)  for  one,  two,  or  three  dimensions, 
reflectively.  Determination  of  the  required  weighting  to  realize  the  ideal 
pattern  (34)  in  different  numbers  of  dimensions  is  taken  up  in  a  later  section 
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If  o  <  -  the  pattern  (28)  has  increasing  si  delobes  as  u  increases;  see 

(30).  Therefore,  a  >  -  are  the  only  cases  of  practical  interest  for  pattern 
(28). 

Plots  of  pattern  (28)  are  given  in  dB  in  figures  8-13  for  a  =  2,  1.5,  1, 
.5,  0,  -.5,  respectively,  for  various  values  of  B.  The  program  is  listed  in 
appendix  D.  The  larger  values  of  a  realize  the  more  rapid  decay  of  sidelobes, 
but,  on  the  other  hand,  have  wider  mainlobes.  Figures  8-13  indicate  the 
necessary  tradeoffs  between  mainlobe  width,  sidelobe  decay,  and 
mainlobe-to-sidelobe  ratio  that  must  be  considered  in  any  weighting  selection. 

A  small  chart  in  the  upper  right  quadrant  of  each  figure  indicates  some 
allowable  values  of  u  and  v  that  apply  to  that  figure.  For  example,  in  figure 
8,  the  pattern  for  a  =  2  applies  to  all  the  following: 

1  3 

u  3  -  (one  dimension)  with  v  3  j 

or 

u  »  0  (two  dimensions)  wi th  v  3  1, 
or 

u  = (three  dimensions)  with  v  =  .  (37) 

When  we  come  to  figure  11,  for  a  =  .5,  however,  the  case  of  u  *  -j,  v  *  -1  has 
an  asterisk  next  to  the  v  3  -1  entry.  The  reason  for  this  is  that  the 
integral  (23)  leading  to  pattern  g(u)  was  convergent  only  for  v  >  -1,  and  now 
we  are  trying  to  violate  that  condition.  A  similar  cautionary  note  is 
indicated  in  figures  12  and  13;  In  fact,  all  three  cases  in  figure  13  violate 
the  condition  v  >  -1.  Oe^jite  this  preclusion,  we  shall  find  later  that  the 
required  weighting  does,  in  fact,  have  the  form  (22)  for  the  corresponding  v 
values  given  in  figures  11  -  13,  but  requires  generalized  functions  with 
singularities  at  the  endpoint  s  3  1  of  the  Interval.  This  extension  to  v  <  -1 
is  desirable  and  Important  because  realization  of  the  ideal  pattern  (figures  7 
and  13)  requires  values  for  v  in  thi  s  region,  regardless  of  the  number  of 
dimension  s. 
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Figure  12.  Pattern  in  dB  for  a  =  0 
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FIRST  NULL  LOCATION 

Let  za  be  the  smallest  nonzero  null  location  of  Ja(z);  i.e. , 

Ja(za)  =0.  (38) 

A  short  list  of  {za”j  is  given  below  in  table  2.  Then,  by  use  of  (17),  the 

Table  2.  First  Zero  of  J^z) 


a  Z 

a 


-0.5 

1.5708 

0.0 

2.4048 

0.5 

3.1416 

1.0 

3.8317 

1.5 

4.4934 

2.0 

5.1356 

2.5 

5.7635 

3.0 

6.3802 

3.5 

6.9879 

4.0 

7.5883 

4.5 

8.1826 

first  null  location  of  pattern  g(u)  in  (28)  is  at  u0,  where 


The  results  in  figure  14  display  the  first  null  location  as  a  function  of 

8,  for  various  values  of  a.  For  large  8,  u„  behaves  as  8  +  -izVB.  By  the 
Identification  of  a  as  »  +  v  +  1,  this  curve  applies  to  any  nunber  of 
dimensions  and  to  whatever  value  of  v  Is  selected  in  weighting  w(  s)  of  (22). 
The  curves  Indicate  that  the  first  null  location  u0  is  monotonic al ly 
Increasing  in  both  B  and  a. 
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LOCATION  OF  PEAK  OF  FIRST  SIOELOBE 


8y  use  of  (28)  and  (21),  we  obtain  derivative 


(40) 


Therefore,  reference  to  (38)  reveals  that  the  location  of  the  peak  of  the 
first  sidelobe  of  g(u)  occurs  vrtiere  g'(up)  =  Q;  i.e. 


(41) 


If  we  employ  more  explicit  notation  in  (3y)  and  (41),  we  can  express  the  first 
peak  location  in  terms  of  the  first  null  location  according  to 

Up( B,a)  *  U0(8,a+1).  (42) 

Thus  all  the  results  in  figure  14  can  be  applied  directly  to  the  first  peak 
location.  For  example,  (42)  yields 

1^,(8,  -  -j)  *  uq(B,  -j);  (43) 

thus  the  third  curve  from  the  bottom  in  figure  14  gives  the  location  of  the 
peak  of  the  first  slde’obe  when  a  *  -  •£. 

PEAK  SIDELOBE  LEVEL 

The  value  of  voltage  pattern  g(u)  at  location  (41)  gives  the  level  of  the 
peak  sidelobe: 


9(v  ^.(j/4**)  ■#.(*« 0- 


(44) 
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Since  the  origin  value  of  the  pattern  is 

9(°)  =  <j^(3),  (46) 

the  voltage  ratio  of  peak  sidelobe  level  to  raainlobe  height  is 


9<u»'  JA* 0 

a 


(46) 


Relation  (46)  is  plotted  in  dB  in  figure  15  as  a  function  of  B,  for  various 
values  of  a;  i.e. , 


dB  =  10 


log 


1 


(B) 


(47) 


The  peak  sidelobe  level  decreases  monotonically  with  increasing  B,  but  has  no 
simple  behavior  versus  a  except  for  very  small  B  or  very  large  B. 


The  results  of  these  last  two  figures  are  combined  in  figure  16,  where  we 
plot  the  peak  sidelobe  level  in  dB  versus  the  first  null  location  uQ,  These 
latter  curves  are  virtually  linear  over  a  wide  range.  If  we  disregard  the 
sidelobe  decay  rate,  the  most  desirable  region  of  this  figure  is  in  the  lower 
left  quadrant,  i.e.,  snail  u0  and  very  negative  dB.  However,  the  closest  we 
can  get  (from  our  family)  Is  the  a  =  -.5  curve,  which  is,  in  fact,  the  ideal 
pattern;  see  (34)  and  figure  7.  Furthermore,  the  sidelobe  decay  rate  is  then 
0  dB /octave.  Higher  sidelobe  decay  rates  are  attained  by  moving  toward  the 
upoer  right  quadrant  of  the  figure;  for  example,  the  a  =  3.5  curve  has  a  6a+3 
*  24  dB/octave  sidelobe  decay  rate.  This  figure  furnishes  a  very  compact 
dl^lay  of  the  Important  interrelationships  that  occur  between  the  fundamental 
features  of  peak  sidelobe  level,  malnlobe  width,  and  sidelobe  decay  rate,  and 
allows  for  a  quick  tradeoff  comparison  of  alternatives. 
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Figure  16.  Peak  Sidelobe  Level  Versus  First  Null  Location  of  Pattern 

9<„)  (fiV) 
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IDEAL  WEIGHTING  AND  PATTERN 

It  was  shown  in  figure  7,  in  the  previous  section,  that  the  pattern  g(u)  for 

a  =  -  takes  on  a  particularly  desirable  behavior,  namely,  a  narrow  mainlobe 
width  and  a  large  ratio  of  mai  nlobe-to-peak-  si  delobe.  However,  figures  11  -  13 
indicated  that  realization  of  some  patterns  was  apparently  not  possible  in 
certain  dimensions  because  we  were  violating  the  condition  on  parameter  v  <n 
weighting  (22)  that  allowed  for  convergence  of  integral  (23).  Here  ve  will 
address  the  more  general  problem  of  how  to  realize  pattern  (23), 

9{u)  |u+v4l  (V-*2)  for  a11  u*  (48) 

for  any  u  >  -1,  but  without  the  current  restriction  of  v  >  -1  in  weighting 
(22).  This  procedure  will  of  course  require  a  different  and  more  general 
weighting  than  (22),  and  will  furnish  solutions  to  the  asterisked  cases  in 
figures  11  -  13. 

The  solution  for  the  required  weightings  to  realize  (43)  for  any  n  >  -1  i s 
conducted  in  appendix  E.  All  the  weightings  are  zero  for  s  >  1,  as  desired; 
their  values  for  0  £  s  S  1  are  listed  below.  From  (E-14), 

for  v  =  -1 ;  ( 49) 


for  -2  <  v  <  -1;  (50) 

] 

and  from  (E-39), 

Ms)  -  ijUs-l)  -  S'U-1*  for  v  =  -2.  (51)  | 

I 

l 

The  extended  range  for  v  <  -21s  given  In  (E-35)  and  (E-36).  Weighting  (49) 
requires  a  generalized  function,  namely,  a  delta  function,  with  its 
singularity  located  at  s  *  1.  Weighting  (51)  requires,  in  addition,  the 
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derivative  of  a  delta  function,  located  at  s  =  1.  The  intermediate  cases  in 
weighting  (50)  require  a  generalized  function  with  its  singularity  located  at 
s  =  1;  interpretation  and  approximation  of  this  generalized  function  is  given 
i n  apoendix  F. 

It  can  be  observed  from  ( 49)  -  ( 51 )  that  the  leading  term  of  w(  s)  is  exactly 
what  would  have  been  obtained  from  initial  weighting  (22)  by  substituting  the 
appropriate  value  of  v;  here  we  are  using  I_n(z)  =  In(z)  15;  9.6.6J. 

However,  the  price  of  crossing  the  “natural  boundary"  at  v  =  -1,  which  was 
originally  required  for  (23),  is  a  generalized  function  with  its  singularity 
located  at  s  -  1.  And  the  further  we  go  below  v  =  -1,  the  more  singular 
becomes  the  required  generalized  function;  these  points  are  elaborated  upon  in 
appendix  F. 

The  explicit  assignment  of  v  values  in  ( 48 ) - ( 51)  leads  to  table  3  for  the 
weighting-pattern  pairs.  With  regard  to  application  of  ( 48 ) -( 51 )  to  the  array 

Table  3.  Weighting- Pattern  Pairs;  u  >  -1 


Weighti  nq _  Pattern 


processing  application  in  various  nunbers  of  dimensions,  we  have  table  4  for 
the  required  weightings,  where  we  have  specialized  the  values  of  n.  In  all 
cases,  the  pattern  realized  is  the  ideal  one  of  (34): 

gj( u)  »^(ju2-B2’)=  ^  cos^uW)  for  all  u.  (52) 
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Table  4.  Required  Weighting  for  Ideal  Pattern 


Number  of  Dimensions 

Value  of  u 

Required  Weighting 

1 

1 

'  7 

(49) 

2 

0 

( 50 )  wi  th  v  =  -  ^ 

3 

1 

7 

(51) 

The  weighting  given  by  (49)  for  one  dimension,  namely,  u  =  -1/2,  has 
already  been  presented  by  van  der  Maas  C93-  However,  the  application  of  (49) 
to  the  realization  of  (48)  for  any  u  >  -1  i s  new.  Additionally,  all  the 
results  in  (50)  and  (51)  for  any  u  >  -1,  and  their  application  to  two-  and 
three-dimensional  array  processing  in  table  4,  are  new.  An  approximation  to 
the  ideal  pattern  in  two  dimensions,  namely,  weighting  (50)  with  v  =  -1.5,  is 
treated  i s  detai  1  in  appendix  F. 
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SOME  OTHER  WEIGHTINGS 

Another  candidate  class  of  weightings  for  consideration  is 


w(  s) 


for  0  <  s  <  I, 


(53) 


of  which  (49)  and  (51)  are  repre  sentati ve  examples,  exclusive  of  the 
generalized  functions.  This  class  is  somewhat  similar  in  fonn  to  the  earlier 
case  in  (22).  Substitution  of  (53)  in  (15)  yi eld s  pattern 


g(u)  =  j'  ds  s(±J  J(us)  B2v  rfv(s^77).  (54) 

Thi  s  integral  converges  at  s  =  0  for  u  >  -1,  but  needs  no  restriction  on  v 
whatsoever. 


To  our  knowledge,  evaluation  of  (54) 

i  s  not 

possible 

in  closed  form 

for 

general  v;  however,  the  following  cases 

are  evaluated  in 

appendix  G: 

g(u) 

2  2,N\ 

1  -u ; 

for 

v  =  0; 

(55) 

g(u) 

=^^/u2-B2’)  -  JfM(u) 

for 

v  =  1; 

(56) 

g(u)  = 

■n 

(u)  for 

v  =  2. 

(57) 

Numerous  jpecia' 

1  cases  for  one,  two,  or 

three  dimensions 

are  available 

from 

(55) -(57)  by  setting  u  *  -  •£,  0,  or  -j,  respectively. 
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As  an  example,  if  we  take  (54)  and  (56)  for  v  s  1,  «  have 


s 


ds  s(-J  J  (us) 


B 


1 


0  fl-  s2 

Addition  of  a  J  function  immediately  yields,  for  u  >  -1 
1  r 


J*f3)  '581 


ds  SW  JU(US) 


B 


.{.P)  »f(  s-l! 


(59) 


which  has  already  been  presented  in  (49)  and  table  3.  A  similar  combination  of 
(54)  and  (57)  yields  the  results  of  (51)  and  the  bottom  entry  of  table  3. 
Additional  results  for  v  =  3,  4,  ...  are  possible  via  the  method  of  appendix  6. 


One  other  two-parameter  family  of  weightings  that  affords  a  closed  form 
pattern  is  available  from  [10;  6.688  1],  by  identifying  u  =  1,  x  =  u,  z  =  iB, 
cos  t  =  s,  and  by  using  [5;  9.6.3  and  9.1.35]: 


This  result  is  restricted  to  the  line  array.  The  weighting  is  continuous  at 
s  =  1  if  v  >  1,  and  the  pattern  (60)  decays  at  3+3 v  dB/octave.  How  good  thi s 
pattern  is  has  not  been  pursued. 


All  the  above  results  have  been  aimed  at  getting  closed  form  results  for 
the  pattern;  however,  this  is  by  no  means  necessary.  One  could  consider  the 
class  of  weightings  (53)  for  any  v,  or  the  class 

exp(-B2s2)  (l-s2)v  for  0  <  s  <  1  (61) 

for  example,  numerically  by  substitution  in  (12)  or  (15)  and  use  of  some 
integration  rule  like  Simpson's.  Once  the  patterns  have  been  numerically 
evaluated  and  plotted  for  a  sufficiently  broad  range  of  values  of  B  and  v, 
good  candidates  can  be  selected  at  will  and  the  corresponding  weighting,  (61) 
or  (53)  for  example,  easily  evaluated.  For  the  line  array,  this  numerical 
approach  is  readily  accomplished  by  use  of  an  FFT. 
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DISCUSSION 

The  ideal  pattern  was  defined  in  (34)  as  (2/w  ^  cos(^u  2-B^),  and  the 
corre  giondi  ng  weightings  were  given  in  { 49 ) -( 51 ) .  Now  in  the  one-dimensional 
application,  u  =  -1/2,  van  der  Maas  [9]  has  indeed  shown,  by  taking  the  limit 
of  a  Dolph-Cheby shev  discrete  array  design,  that  there  is  no  pattern  with  a 
narrower  mainlobe  for  a  specified  sidelobe  level  (and  vice  versa)  than  (34). 
However,  stn'ctly  peaking,  we  have  not  proven  that  same  result  for  the  other 
values  of  u,  i.e. ,  other  numbers  of  dimensions.  Instead,  we  have  adopted  (34) 
as  an  ideal  pattern  and  shown  that  it  can  be  realized  by  finite-support 
weighting  functions  with  a  generalized  function  whose  singularity  i  s  centered 
at  the  edge  of  the  array.  Conceivably,  there  might  be  a  different  weighting 
that  would  realize  a  pattern  that  gets  further  into  the  left-lower  quadrant  of 
figure  16.  However,  we  conjecture  that  this  is  not  possible  and  that  the 
leftmost  curve  in  figure  16  is  the  ultimate  attainable  region  for  any 
weighting  in  any  ntmber  of  dimensions. 

SUMMARY 

We  have  presented  a  two-parameter  class  of  Bessel  weighting  functions  that 
have  a  closed  forw  pattern  with  controllable  mainlobe  width,  mainlobe-to- 
peak-sidelobe  ratio,  and  sidelobe  decay  rate.  These  results  have  application 
to  arrays  in  one,  two,  or  three  dimensions.  In  addition,  the  ideal  pattern 
and  the  corregjondi ng  weightings  required  in  various  numbers  of  dimensions 
have  been  derived  and  presented.  Where  a  generalized  function  is  required,  a 
method  of  approximating  it  has  been  presented  and  illustrated  by  a  numerical 
example.  Various  weightings  already  extant  in  the  literature  were  shown  to  oe 
special  cases  or  limiting  cases  of  the  general  results  given  here. 
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Appendix  A 

DERIVATION  OF  RESPONSE  OF  PLANAR  ARRAY 

Let  the  receiving  array  lie  in  the  x-y  plane,  and  let  a  plane  wave  of 

wavelength  x  arrive  at  polar  angle  da  and  azimuthal  angle  ea;  see  figure 
A-l.  Then  if  the  time  of  arrival  of  the  plane  wave  at  the  origin  is  denoted 


z 


Figure  A-l.  Geometry  for  Planar  Array 

as  0,  the  time  of  arrival  at  a  general  point  x,y  in  the  plane  of  the  array  is 


x  cos  e,  +  y  sin  e. 

ra  =  -sfn<Sa - ; - 


(A-l) 


where  x  =  c/fQ>  c  fS  the  speed  of  propagation,  and  f0  is  the  frequency  of 
the  plane-wave  arrival. 


To  look  in  direction  ,  the  receiving  array  should  employ 

steeri  ng-delay 


Xx  -  -sin  4  x  cos  a,  +  y  sin  e, 


(A-2) 
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at  the  general  point  x,y.  If,  also,  weighting  w2(x>y)  1S  used  in  the 
receiving  array,  the  array  output  voltage  response  at  frequency  f  j  s 


g  *  [fdx  dy  w2(x,y)  exp  [ji 2»f o(Ta+ 

=  JJdx  dy  w2(x,y)  exp Jj-i 2w  j  P1  -  i2*  ^  P2J  ,  (A-3) 


where  we  define  angle  functions 

P1  =  sin$cos^"s"in!(bcos^i »  P2  =  sin^  sin^ -sini6asinea,  (A-4) 

and  where  the  integration  is  carried  out  over  all  x,y  where  weighting  w2  /  o. 
Thus  the  planar  array  can  have  arbitrary  geometry  in  the  x,y  plane;  equations 
(A-3)  and  (A-4)  are  general  results  for  the  array  response. 

If  weighting  w?  contains  an  impulse  at  x0,  y0,  then  we  have 


w2(x,y)  =  $  (x-x0)|(y-y0). 


(A-5) 


with  array  response 


g  =  exp  jji  2*  ^  pj-l  2tt  ^  P2j  , 
which  never  decays  In  amplitude  with  increasing  angle. 
SPEC IALIZATION  TO  CIRCULAR  ARRAY 


( A-6) 


As  a  special  case  of  the  above,  consider  a  pi anar-circular  array  of  radius 
R  with  weighting  Independent  of  angle;  i.e. , 


for  x  +y  <  R 
otherwt  se  j 


(A-7) 
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Then  regionse  (A-3)  becomes 
R 

9  ■  Idr 
0 


*  r 

(  de  w^(  r)  exp  ^-12»y^P 

1  cose  +  P2  sine) 

=  2,  J  dr  r  wjtr)  ^l2+p?)* 


And  from  (A-4),  we  have 


2  2*  T  2  2  ~)^ 

VP1  +P2  =  S1'n  ^  +  sin  **a  "  2sin^  si  n^  cos  (^-ea)| 


le. 


sin^e  -  sindae 


Mow  let 

r  «  Rs 

in  ( A— 8) .  There  follows,  for  re^onse  (A-8), 

1 

g(u)  »  J  ds  s  J(us)  »i(s), 
0 

where 


li!  2l 


si  e  -  si  n<Jfle 


and 


w(  s)  *  2»R^Wj(Rs)  for  0  <  s  <  1. 


(A-8) 


(A-9) 


(A-ld) 


{ A-ll) 


( A-l  2) 


(A-U) 
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Thus  the  voltage  response  g(u)  is  given  by  a  zero-th  order  Bessel 
transform  over  (0,1)  of  normalized  weighting  w(s).  Dimensionless  parameter  u 
incorporates  the  received  wavelength  x,  the  array  radius  R,  and  the  various 
look  and  arrival  angles. 


40 
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Appendix  B 

DERIVATION  OF  RESPONSE  OF  VOLUMETRIC  ARRAY 

Consider  a  plane  wave  of  wavelength  x  arriving  at  angles  da>  ea;  see 
figure  A-l.  The  time  of  arrival  (relative  to  the  origin)  at  a  general  point 
x,y,z  in  the  volumetric  array  is  [11;  eq.  3J 


(8-1 ) 


T  3  -  ^fx  cose  sinri  +  y  sine,  sin  d  +  z  cos/6  1. 
acLaa  a  a  aj 

To  look  in  direction  *1,  ,  e^ ,  the  delay  used  at  location  x,y,z  should  be 


^  =  c  [x  Cos^  s1n^  +  y  sin^  sin  ^  +  z  cos^J.  (B-2) 

The  response  of  a  weighted  array,  at  frequency  fQ>  is  then 

g  =  £^dx  dy  dz  w3(x,y,z)  exp  [-i 2irfo(Ta+Tj)] 

=  dx  dy  dz  W3 (x,y,z)  exp  i^(x P^y ?z+z ,  (8-3) 


where  W3  i  s  the  weighting  and 


Pj  =  cos  sin  ^  -  cos  ea  sin  6a, 
P?  =  sin  ^  sin  ^  -  sin  eg  sin  da, 


?3  =  cos  ^  -  cos  ^a< 


(8-4) 


The  Integration  In  (3-3)  is  over  all  x,y,z  where  W3  f  0.  (b-3)  is  the 

general  result  for  any  time-delay  steered  volumetric  array. 
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SPECIALIZATION  TO  SPHERICAL  ARRAY 

Let  veighting  W3  depend  only  on  the  distance  from  the  array  center;  that 
i  s. 


v^U.y  ,z)  =  Wj 


2  2  21  222  2 
1  +y  +z  )  for  x  +y  +z  <  R  , 


(3-5) 


where  R  is  the  sphere  radius.  Then  the  voltage  response  of  the  array  is,  from 
(8-3)  and  (B-5), 


1 

1 


(B-7) 


IK  b/61 


where  we  define 

0. 

Then  the  response  is,  upon  substitution  of  (B-7)  in  (8-6), 


(B-8) 


I 


g  =  4*  dr  r  w,  ( r) 


sin(2"lQ) 


2^0 


{ B-9 ) 


Now  let  r  =  Rs;  then  the  response  can  be  written  as 


1  . 

& 

g(u)  =  J  ds 

2)s  Msl, 

it/  u 

(B-10) 

0 


where 

w(s)5  ( 2^)3/2  R3w1(Rs)  ,  (B-li) 

and  dimensionless  parameter 

u  s  2i£q.  (B-12) 

The  quantity  Q,  involving  the  look  and  arrival  angles,  can  be  expressed  as 


Q  *  J^2-2  cos^  cos«ia-2sin^  sin»*a  cos(o^  -eajj.  (3-13) 

Thus  the  general  result,  (B-10),  for  the  volumetric  array  is  given  by  a  sine 
transform  of  the  normalized  weighting  w. 


i 


I 
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Appendix  C 

MONOTONICITY  OF  BESSEL  WEIGHTING 
We  investigate  here  the  behavior  of  weighting  (22). 


w(  s) 


Let  y  ■  3^1-  s^;  then  (C-l)  becomes 


for  0  < 


w(s)  =  B‘2v  yv  lv(y). 


Now  [5;  9.6.28] 


which  is  positive  for  v  >  0,  y  >  0;  see,  for  example,  15; 
i  s  monotonic  ally  increasing  in  y  if  v  1  0;  therefore  (C-l 
decreasing  in  s  if  v  i  0. 


TR  o761 


s  <  1.  (C-l) 


(C-2) 


( C— 3 ) 

9.6. lOj.  Tnus  (C-2) 
i s  monotonically 
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Appendix  D 

PROGRAM  FOR  CALCULATION  OF  PATTERN  (28) 


10 

ft  1  ph  a=-  .  5 

20 

PLOTTER  IS  "98720" 

30 

LIMIT  25, 175. 35,245 

40 

OUTPUT  705; " V S 5 " 

50 

SCALE  0. 6*-PI , -70, 0 

60 

GRID  PI ,  10 

70 

PE  HUP 

30 

FOR  B  =  2  TO  8  STEP  2 

90 

B2=B*6 

1  00 

F  =  FN I rig  nu<  A  1  pha,  B  > 

1  10 

FOR  IJ  —  0  TO  B  STEF'  .05 

120 

'.'=FN I nu  tiu 1  A1  pha,  SQR<B 

2-U*U 

1  30 

PLOT  U,  20*LGT-;ABSCY/F> 

> 

1  40 

NEXT  U 

150 

FOR  U=B  TO  6*P I  STEP  . 

05 

160 

V = F  N  J  n  u  x  n  u  <  A  1  p  h  a ,  S  Q  R 1  U  *  U  -  B  2 

1  70 

PLOT  U ,  2 0 - L G T  <  ABS  <  Y.  F  > 

> 

180 

NEXT  U 

190 

PE  HUP 

2O0 

NEXT  B 

213  EHD 
220  1 

230  DEF  FNG  arnma1  X  >  1  GAMMA<X>  Mia  HART,  page  279 ,  #'5231 
240  N= I  NT  <  X > 

250  R=y,~N 

260  IF  <  N , 0  ‘  OR  <R< >»>  THEN  290 

270  PRINT  "FNGammaOC  IS  NOT  DEFINED  FOR  X  =  " ;  X 
2S'0  STOP 

290  IF  R  0  THEN  320 
300  G  am  m  a2 = 1 

310  GOTO  260 

320  P  =  3 . 36954359131 +R*f' 1 . 09850630453+R* " . 1 42928007949+R*3 . 930 1 3464 1 8SE-2 > 

330  P  =  4 3 . 94  1  0209 1 89  +  P*  <  22 . 96308008 36+R* 12.8021696 1  12^R*P  ;•  . 

340  0=43. 94 1 0209 191 +R*  <  4. 39050474596-R# < 7. 15075063299-R) ) 

350  Gamma2=PG  !  GAMMA <2+R>  for  O  R  1 

360  IF  N  2  THEN  400 

370  IF  N< 2  THEN  450 

380  Gamma3 Gam ma2 

390  GOTO  500 

4u0  Gamma=Gamma2 

410  FOR  K= 1  TO  N-2 

420  Gamma*G  amnia*  <  X-K 

430  NEXT  K 

440  GOTO  500 

450  R* 1 

468  FOR  K=0  TO  1-N 
470  R = R * < X  +  K.  ) 

480  NEXT  K 

490  G  am  ma= G  am  m  a2 'R 

500  RETURN  Gamma 

510  FNEND 

520  ! 


47 


TR  6761 


5  30  BEF  FN Jnuxnu  <  Nu , X  >  !  JriuOi.:'  "i‘'U 

540  IF  hBS<  ::  >  1  THEN  65U 

550  ft*. 797834560303 

560  IF  Mu  =0  THEM  RETURN  FNJoOO 

570  IF  Nu«.5  THEM  RETURN  ft*S  I N  X 

560  IF  Hu  =  - .  5  THEM  RETURN  ft*COSOO 

590  IF  Mu- 1  THEM  RETURN  FMJ 1  iX  1  .-'X 

600  IF  Mu  =  - 1  THEN  RETURN  -FM  J  1 '  X  >  *X 

610  IF  Mu*  1.5  THEN  RETURN  ft  *  <  3 1 N '  X  •  -  X  *  C  0  6  :  •  X  3 

610  IF  Hu— 1.5  THEN  RETURN  -  ft  *  <  X  *  S I  N  ■,  X  >  +  C  U  3 1  X  1  ■ 

6  30  IF  Mu  *2  THEN  RETURN  ( 2*FN J 1 < X > -X*FN Jo ■  •  : ■  3 

640  IF  Nu*-2  THEN  RETURN  <2*FN J 1 < X >-X*FN Jc •  : '  > 

650  ft  =  N 

660  IF  ■  I N  T 1  fl  > <  H  1  OF:  <ft>  =  0.)  THEN  630 
670  K  =(i  =  -Nu 

6 8 0  3  *  T  =  1  1 2  fi*FNG am m a < 8  + 1 ) > 

690  R  =  - .  25*X*X 
700  B  i  9=HBS  <  3  > 

710  FOR  M* 1  TO  10O 
720  T  =  T  *R/  (.  N*'.  N  +  ft  >  > 

7  3  0  '3  =  S  +  T 

740  B i  g=MRX <,  B i  g ,  hBS 6 >  > 

750  IF  flBS<TX«lE-11^0BS<S>  THEN  790 
760  NEXT  N 

770  PRINT  "100  TERMS  IN  FNJriuxnu'.  Nu.  7!  '  FIT  "SHuj:-: 

7 SO  PftUSE 

790  D=  1  2-LGT  (HBS  <.  B  i  g-  S  '  '  NO.  OF  SI  ON  IF  DIGITS 

SOO  IF  K  > 0  THEN  S  *  S  ‘  4  *  R  >  •'  K 


310  RETURN  S 
320  FNEND 

3  30  1 

840  I)EF  FNJoJX ■  1  Jo<X>  ui  a  9.4.1  ”,  9.4.3  of  ftr'S  55 

350  Y*FiB$<X> 

860  IF  V 3  THEN  310 
370  T  *  Y  * Y -  9 

380  Jo*. 0444479-T*' . 0039444-T*. 00021 > 

390  Jo*l-T*'  2. 249999 <'-T4(  1 , 2656203- T*‘  .  3 1 6  :  366- T*  Jo  ■  . 

900  GOTO  970 
9l0  T  *  3  '  Y 

920  Jo*9. 512E-5-T*' . U0137237-T*< . 00072305-T  - . 000 1 4476  • 

9  30  Jo*.  79788456-T*.  7. 7E-7*T*<. .  005 52 740+ T *  Jo 
340  3*. 00262573-T* •  . 00054 1 25*  T*  r . 00029333 - f«.800t  3558  • 

950  T-Y-.  73539316-T*.  .  04  1 66397  +  T*  ■;  3 . 954E-5- T  -  S  •  > 

960  Jo* Jo*CQS< T  •  SQR'  YJ 
970  RETURN  Jo 
980  FNEND 
990  1 

1000  DEF  FUJI  >'  X  '  1  Jl-X,  via  9.4.4  *  9.4.6  ft  MS  55 

1O10  Y*ftB3‘ X  ■ 

1020  IF  <  3  THEN  1070 
1030  T*Y*Y  9 

1O40  J1 *. 0044331 9- T  +  <  . OO031 761 -T*. 0008 1 1 09  ■ 

1050  J  l  «.<*•' .  5-T*'  .  562‘ 9985-T*.  .  2  1 89 3573- T  •» •  .  8  3.>5428'=‘-f* 
1068  GOTO  1130 
1O70  T*  3  i 

1O80  J 1  *  .  00017185-T *<. .  O02495  1  1  -  T  ■  .  O0 1 1  365  3-T* .  88820833 
1098  J  1  *.  79788456+T* ■' J  .  56E-6  +  T* < .  01  659667  +  T*  J  1  •  1 
1100  S*. O06 37379-T ♦ <  . U007434S  +  T<^ . 00079824-T *.800291 66  • 

1110  T * Y - 2 .  356  1 9449*T ■>■<  .  1 24996  1  2  +  T*  <  5 . 65E -5  - T* S  ■  • 

1  120  J  1  *  3GN '  X>*J1*C0S<  T  >  SQR  <  Y  j 
1130  RETURN  J1 
1140  FNEND 
1150  1 
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l 

1  170 
1  1  50 
1  1  90 
1  100 
1510 
1520 
12  50 
1540 
1550 
12*0 
12  70 
1250 
1290 
1  500 
1510 
1520 
1330 
1340 
1250 
13*0 
1370 
1330 
1390 
1400 
14  10 
1420 
1430 
1440 
1450 
1460 
14  70 
1430 
14  90 
1500 
1510 
1520 
1530 
1540 
1550 
15*0 
1570 
1580 
1590 
1*00 
1*10 
1*20 
1*30 
1*40 
1*50 
1  **0 
1*70 
1*30 
1*90 
1  7O0 
1710 
1720 
1730 
1  740 


DEF  Ftlli'iij.  nu'  Hu,  !■!.)  !  I  ru.  ■  nu 

IF  AB3 1  :: 1  THEN  1290 
H- . 39394228040 1 

e=e;:po.  • 

IF  Nu  =  o  .THEN  RETURN  FNIoCX; 

IF  Nu=.5  THEN  RETURN  m*',E-1  E'  X 
IF  N  u  =  -  .  5  THEN  RETURN  A*<E+1  E  • 

IF  Nu  =  1  THEN  RETURN  FN  I  U X  •  /  II 
IF  Nu= - 1  THEN  RETURN  FNI1 

I F  Mu*  1 . 5  THEN  RETURN  (  X- t)4E+«,  :>  1  E  ■  •  X-  3 
IF  Mu=-1.5  THEN  RETURN  ft*<<X-l  >*E-0-.+  .  •  'E  • 

IF  Nu«2  THEN  RETURN  <X*FNIg<X > -2*FNI  1  •  .  3 

IF  Nu=-2  THEN  RETURN  c.X*FNIo<X  ••-2*FHI  1  X> 

M  =  Nu 

IF  ■  INT'  H  ■  H  ■  OR  •  (H,  =0>  THEN  1350 
K  =  A  =  -  H  u 

5  *  T  =  1  /  '  "  A  +  FNG  arnfii  a*.  Fl+  1  >  ’> 

R*. 2  5  *  X  *  X 
B  i  g=AB3  <.  S  > 

FOR  N= 1  TO  100 
T  =  T4R,"CN4(N  +  H>  > 

S  =  3  +  T 

B i  g=MAX ■'  Big, ABS<$>  > 

IF  ABS<T><=1E-1  1 +ABS  <  3  )  THEN  1430 
NEXT  N 

PRINT  "100  TERMS  IN  FN  I  nuxnu  <  Nu ,  X  ■  FIT  "  ;  Mu  5  X 
PAUSE 

D=  1 2-LGT  < AB3  C B i g- 3  > )  !  NO.  OF  SIGN  IF.  DIGITS 

IF  K  > 0  THEN  3  =  3*  < 4*R>-K 
RETURN  3 
FNEHD 
! 

DEF  FN  I  o  < X >  !  Io <. X  )  m U  9 . 3  .  1  «  9 .3.2 

Y*ABS<X> 

IF  Y >3 . 75  THEN  1550 
T=Y*Y  '  14. 0*25 

I o= . 2659732+T *  <  . 03*0763* T * . 00453 1 3  > 

I  o=  1  +  T *  <  3 . 5  1 5322 9  +  T * (3 . 039942 4  +  T * 1  1 . 2u*7 492  +  T -►  I  :  >’>  • 

GOTO  1590 
T  =  3. 75,-Y 

I  o  = .  0  0  9 1  *  2  3  1  —  T  *  ■  ,  02057706— T*1  .  02635537  —  T  *  '■  .  01647*33— T*.  O  O  3  9  2  3  7  7  '> )  > 

I  o=.  39394228 +  T*(  .  01  328592 +  T*<  .  00225  3 1  9-T* •. .  @0 1 575*5- T* I o  >  >  > 

I  o  =  I  o  *  E  X  P  *  Y  >  3  Q  R  <  Y  > 

RETURN  Io 
FNEND 

I 

DEF  FN  II  <  X  >  'll  ■  X  >  w  i  a  9.3.3  8.  9 .3.4 

Y  =  ABS  < X  t 

IF  Y >3.75  THEN  1*90 
T  =  Y+Y  ■'  1  4. 0625 

II  =  . 02*5873  3 +  T*>'  .  0030 1 532 +T* . O80324 1 1  ■ 

I  1=X*<  .  5*T*<  .  37390594+T+C  .  51498869+T*-  .  15034934+T  *  I  1  ■  ,  > 

GOTO  1730 

T  =  3 . 75> Y 

I I  = .  0  1  0  3  1 555- T  +  <: .  02232967-T*'-  .  02895  3  1 2-T*  (  .  y  1  7 37654-T*  .  @0420059  :> > 
II*.  39894223-T  *  >.  .  03988024+T*'  .  O036201S-T*-  .  001*38@1-T*I  1  '  -■ 

I  i*sgn«  ::>*i  i *EXP‘  y:v.  sqr<y> 

RETURN  II 

FNEND 
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Appendix  E 

DERIVATION  OF  WEIGHTING  FOR  IDEAL  PATTER N 
We  want  to  realize  pattern  (23), 

9^)  forallu’  (E-l) 

for  any  u  >  -1,  but  without  the  restriction  v  >  -1,  which  was  required  for 
convergence  of  integral  (22).  $  a(  z)  was  defined  in  (17)  et  seq. 

We  begin  the  derivation  for  the  requi red  weighting  to  realize  (E-l)  by 
substituting  (E-l)  in  (16): 


00 

w(  s)  =  ■ 

0 

1  d""(l)  l( 

for  s  >  0. 

(E-2) 

It  is  important  to  observe  that  we  must  allow  all  s  >  0  in  (E-2);  hopefully, 
when  we  evaluate  w(  s)  from  (E-2),  it  will  be  zero  for  s  >  1. 


Now  we  already  know  from  (22)  and  (23)  that  (E-2)  yields 


w(  s) 


for  0  <  s  <  i 

for  1  <  s 


if  u  >  -1,  v  >  -1.  (E-3) 


Letting  a  =  u+v+1  in  (E-2)  and  (E-3),  and  eliminating  v,  we  have  the  useful 
integral  identity 


For  convergence  of  this  integral  at  u  =  0,  we  require  >i  >  -1,  whereas  for 
convergence  at  u  =<»,  we  must  have  o  >  u;  i.e. ,  -1  <  u  <  a. 
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Mow  we  have  the  relation  [5;  9.1.30J 


u(l)  JM(SU)  =  -  1  Hl{u(lT  Ju-1(su)]  for 


s  >  0. 


(t-5) 


Thus  (E-2)  can  be  exoressed  as 


Wis)  =  -  |£iu  u(h)  J.1(su)^ii4v+1(yuZ-BZ,|  for  s  >  0.  (E-6) 


Appeal  to  { E -4 )  now  reveals  that 


w(  s)  = 


ij.fi M2 

s  ds  )  V  3 


for  0  <  s  <  1 


0  for  1  <  s 


(E-7) 


provided  that  u  >  0,  v  >  -2. 

We  have  succeeded  in  extending  the  range  of  v  from  v  >  -1  to  v  >  -2,  as 
desired,  but  have  apparently  altered  and  restricted  the  range  of  u  from  »  >  -1 
to  u  >  0.  However,  thi  s  last  restriction  is  due  solely  to  the  method  of 
derivation,  and  may  be  restored  to  u  >  -1,  by  observing  that  the  right-hand 
side  of  (E-7)  is  analytic  in  u  (in  fact,  constant),  and  that  the  function  w(  s) 
defined  by  (E-2)  is  analytic*  in  w  for  »  >  -1.  Thus  (E-7)  gives  the  required 
weighting  to  realize  pattern  (E-l),  provided  that 


u  >  -1,  v  >  -2. 


(E-8) 


Hovever,  care  must  be  taken,  in  the  evaluation  of  the  derivative  in  (E-7),  to 
account  for  any  generalized  functions  that  may  be  generated. 


We  are  using  the  fact  that  $a(z)  i  s  an  entire  function  of  a,  regardless 
of  the  value  of  z;  see  (17)  et  seq. 
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In  (E-7),  define  the  function 


Then  there  follows 

w(  s)  =  -  ^  Wj(  s)  for  alls,  if  u  >  -1,  v  >  -2.  ( t-lo ) 

We  now  break  the  region  v  >  -2  into  the  three  subcases  ( i )  v  >  -1, 

{ i i )  v  =  -1,  (iii)  -2  <  v  <  -1. 

( i )  v  >  -1 


As  s-*  1- ,  Wj(  s)  in  (E-9)  approaches  0,  since  v+1  >  0.  Thus  w^I  s)  in 
(E-9)  is  continuous  for  all  s,  and  we  find,  by  the  use  of  [5;  9.6.28J  in 
(E-7),  that 


w(  s)  = 


which  checks  (22),  as  it  must  of  course,  for  v  >  -1. 


(E-ll) 


We  observe  that 


»4  s) 


(l-s)v 

p(  V+1) 


as  S-+1-. 


(£-12) 


Thus  for  -1  <  v  <  0,  there  is  an  integrable  singularity  in  w(  s)  at  s  =  1.  For 
w  >  0,  w(  s)  is  continuous  at  s  =  1,  while  for  v  =0,  w(  s)  has  a  discontinuous 
step  of  value  -1  at  s  *  1. 
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An  alternative  way  of  obtaining  the  ""suit  ( E - 14 )  is  by  taking  the  limit 
of  (£-11)  as  v-*-l.  First  we  observe  that 


since  Iv(z)  iS  an  entire  function  of  v  wien  z  /  0  L 5 ;  9.6.1J.  We  then 
define  a  difference  or  remainder  function  (for  v  >  -1)  as 
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Therefore,  combining  (E-17)  and  (E-21),  we  obtain 


in  agreement  with  (E-14). 

An  alternative  and  simpler  equivalent  form  of  C E - 22)  is 


{ E  -22) 


JJ-i  (V5)Vlv(V?)  -  +  2*(  s>  for  s  *  °*  <£-23> 

The  derivation  of  (E-23)  is  similar  to  that  given  above  in  (E-13)-(£-21) . 

(ii i)  -2  <  v  <  -1 


We  return  to  (£-7)  and  { E-9) .  Observe  that  [5;  9.6. 7] 


s)~ 


MT1 

2”*1r(v+2i 


as  s 


Si  nee  we  now  have 


-1  <  v+1  <  0, 


(£-24) 


(E-25) 


there  is  an  integrable  singularity  in  w^{  s)  at  s  =  1.  Thus  the  derivative 
in  { E - 10 )  will  generate  a  generalized  function  with  a  singularity  located  at 
s  *  1.  We  handle  this  case  by  defining  an  auxiliary  function 


'M 


v+l 


A{  s)  s 


2v+lp(  v+2) 


for  0 


<  s  <  l) 


for  1  <  s 


(£-26) 
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Then  using  (£-9),  the  difference 


for  0  <  s  <  1> 


wj(  s)  -  A(  s)  = 


is  a  continuous  function  of  s;  in  fact, 


r(v+2) 

0  for  1  <  s 


|(E-27) 


J 


B 


'AzA 


v+2 


W,  (  s)  -  A{  s)~  — , 

1  2V  3P(v+3) 


as  s-*l-. 


(E-28) 


which  approaches  zero  since  v  >  -2.  Therefore,  the  required  weighting  w(  s)  in 
(E-10)  can  be  expressed  as 


Ms)  =  -  ~  [wl(  s)  -  A(  s)  +  A(  S)} 

=  D(  s)  +  G(  s) ,  (E-29) 

where  both  D( s)  and  G( s)  are  zero  for  1  <  s.  The  difference  function  D( s) 
possesses  an  integrable  singularity  at  s  =  1;  In  fact  (recall  ( E-25) ) , 

q2 A  .21 V+1 

D(s)^  - ' -  as  s-»l-.  (E-30) 

2  v*2nv+2) 

The  last  tenn  in  w(  s)  in  (E-29)  is  a  generalized  function;  from  (E-2b), 
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‘  TUT  A(s) 


2vrT(v+2) 


2vP{v+l) 


for  s  >  0, 


where  the  sub  6  denotes  a  generalized  function.  Here  step  function 


U(t) 


0  for  t  <  0 


1  for  t  >  0 


Combining  (£-29)  and  (£-31),  the  required  weighting  is 


(£-31) 


(£-32) 


iv  8  /  A  '  2vr(v+D  2vr<vs 


for  0  <  s  <  1 


0  for  1  <  s _ J 

(£-33) 


Me  have  now  completed  the  consideration  of  the  three  subcases  delineated 
under  (E-10).  We  now  wish  to  extend  v  to  values  that  are  equal  to  and  less 
than  -2,  so  that  we  can  handle  the  volumetric  array  discussed  in  (37)  et  seq. 


Me  return  to  (£-6)  and  employ  (£-5)  again: 


*  S)  *  '  jf dU  ["  111  u(l)  Ju-2(sujj]j 

*  -\A  j"'ld!  J  duu(l)  V2(su)jLv+l  (jfi 


U  +  v+l\ 


LW' 


u  - B  )}}  for  s  >  0. 


(£-34) 
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Recourse  to  { E -4)  yields 


forQ  < s  <  i 

O  for  1  <  s  I 


provided  that  v  >  -3,  u  >  1.  However,  the  last  restriction  on  u  may  be 
restored  to  u  >  -1,  by  the  argument  under  (E-7).  Thus  (£-35)  gives  the 
required  weighting  to  realize  pattern  (E-l),  provided  that 


(£-35) 


u  >  -1,  v  >  -3. 


(E-3o) 


The  only  subcase  of  (E-35)  that  we  consider  in  detail  is: 


(iv)  v  -  -2 


We  now  can  write  (E-35)  as 


w(s)  -  -  4^|l0(Blf7]u(l-s)  'j. 


Then  [5;  9.6.27] 


!  -  ~  |(l-s2)'i(-2s)\](l-s) 

— ; ^  I.  (BV^^JUU-S)  +S(S-1). 


iio(0)[-S(s-ii] 


(E-37) 


(E-38) 


Therefore  [5;  9.6.28], 


59 


TR  6761 


60 


TK  o761 

Appendix  F 

APPROXIMATION  OF  A  GENERALIZED  FUNCTION 

Equation  (50)  or  (E-33)  gives  the  required  weighting,  for  -2  <  v  <  -1,  to 
realize  pattern  (48)  In  terms  of  a  generalized  function  which  is  difficult  to 
interpret.  Here  we  address  thi  s  i ntrepretation  by  means  of  an  approximation 
to  the  generalized  function  (E-31).  We  begin  by  approximating  the  (singular) 
auxiliary  function  A(s)  defined  in  (£-26)  by  an  ordinary  function  Ag(  s)  ^  an(j 
then  derive  an  approximation  to  generalized  function  G( s)  of  (E-31)  according 
to  the  same  rule,  namely. 


G  (s)  = 

e 


1  d  .  .  . 

-  7  -77  A  {  s) 
S  dS  e 


( F-l  i 


In  particular,  consider  Ae(s)  as  shown  in  figure  F-l;  that  is,  Ae(  s) 
is  still  given  by  ( E - 26 )  for  0  i  si  1-e,  but  then  tapers  linearly  to  zero  at 
s  =  1  in  order  to  be  everywhere  continuous.  The  height  H  of  A£( s)  at  s  =  1-e 
i  s,  from  ( E -26) , 


_  e V+1( 1-e/ 2) U+1 


(>  0). 


(F-2) 


For  snail  e,  we  have 


rfv+2)  as  e**0+’ 


which  tends  to  +oo  as  e-*0+,  since  we  have,  from  (E-25), 


-1  <  v+1  <  0. 


(F-3) 


( F  —4 ) 


The  result  of  applying  (F-l)  to  figure  F-l  is  shov*i  in  figure  F-2.  The 
large  positive  pulse  in  (1-e,  1)  has  height  proportional  to 


H  .  eV(  l-e/2)v41 


nw-2) 


as  e-*0+. 


( F— 5 ) 
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Figure  F-2.  Approximation  to  Generalized  Function  G(s) 
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which  is  tending  to  infinity  since  -2  <  u  <  -1.  The  area  of  this  positive 
pul  se  i  s 


P  U  U  _  V  ^ 

C  =  \  ds  —  =  -  -  1  n{  1-e) ~  -^7 — as  e-*0+.  ( F-6 ) 

J  eS  e  P(v+2) 

1-e 

Thi  s  area  i  s  also  tending  to  +«®  as  e-*0+;  recall  ( F -4 ) .  The  rate  of  increase 

of  C  is  greater,  the  closer  v  is  to  -2.  {For  v-»  -1,  H-*  1,  C-*  ln(l-e);  thus 
area  Ol  as  e-»0+.  This  is  the  unit  area  impulse  presented  in  (49)  for 

v  =  -1.) 


fl 


Figure  F-2  is  one  approximation  to  generalized  function  ii(  s)  defined  in 
(E-31).  Its  most  important  feature  is  the  impulsive-like  positive  pulse  near 
s  =  1.  An  alternative  approximation  is  afforded  in  figure  F-3,  where  the  area 
C  of  the  impulse  at  s  =  1  is  given  by  (recall  (F-6)) 

„  v+1 

The  notation  used  in  (E-31)  for  the  generalized  function. 


G(  s) 


JtK 

2vr(v-t-i) 


(f-8) 


conceals  the  positive  impulsive  behavior  at  s  =  1  that  the  series  of 
approximations  in  figures  F-l  through  F-3  indicate  must  be  present.  In  fact, 
(F-8)  is  negative  for  0  <  s  <  1,  by  reference  to  (F-4). 

The  alternative  approximation  we  obtain  to  weighting  (50)  is  then,  from 
(E-29)  and  figure  F-3, 
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Figure  F-3.  An  Alternative  Approximation  to  Generalized  Function  G($) 


Figure  F-4.  An  Approximation  to  Weighting  (50) 
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we(  s)  =  0{  s)  +  Ge(  s) 


The  plot  in  figure  F-4  illustrates  this  approximation.  The  “transition 
oortion"  in  (1-e,  1),  which  is  the  bottom  line  of  (F-9),  is  singular  at  s  =  1- ; 
however,  this  is  an  integrable  singularity,  as  may  be  seen  by  reference  to 
{ E -29)  and  (£-30).  The  impulse  at  s  =  1  i s  of  finite  area  C  given  by  (F-7). 

As  e  — »0+,  area  C  of  the  impulse  tends  to  infinity;  see  (F-7)  and  (F-4). 
However,  ttie  area  under  the  main  portion  of  the  approximation  w_(s) 
precisely  cancel  s  this  singular  behavior;  that  is,  as  e-*0+. 


which  is  -C.  Since  the  c..*ea  under  the  transition  portion  is 


(F-10) 


d  s  D(  s) 


B2(1-s)v+1 


bV+2 


as  e-*0+, 


( F-ll) 


which  tends  to  0  as  e-»0+,  the  area  under  approximation  wf(s)  remains  finite 
as  e-*0+.  Indeed  it  must  remain  finite,  because  Hankel  transform  (12)  or  (15) 
must  remain  finite  in  order  to  realize  pattern  (48),  which  is  entire  in  u,  v, 
u,  and  6. 
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Thus  our  final  approximation  to  w(  s)  is  simply 


(F-  12) 


and  is  shown  in  figure  F-5.  It  is  very  important  to  observe  that  the  simple 
expedient  of  approximating  w(  s)  by 


( F-13) 


is  totally  inadequate  because,  as  e-*0+,  the  area  under  tne  cusp  at  s  =  1-e 
tends  to  infinity,  and  cannot  possibly  yield  an  entire  function  for  the 
pattern.  The  impulse  i  s  necessary  to  compensate  for  the  singular  behavior  of 
(F-13)  near  s  =  1;  it  allows  us  to  realize  the  "finite  part"  of  the  Hankel 
transform  of  (F-13)  for  e  =  0. 


(As  v-*-l,  the  value  of  C  in  (F-7)  tends  to  the  finite  value  1,  which  is 
the  impulse  in  (49).  And  as  v-»-2,  the  doublet  of  (51)  could  probably  be 
extracted  as  a  limit  from  (F-12);  this  procedure  has  not  been  pursued.) 


The  result  of  using  approximation  (F-12)  with  (F-7),  for  v  =  -1.5  and 
B  =4,  is  displayed  as  the  patterns  in  figures  F-6  through  F-8  for  e  =  .1, 

.01,  .001,  reqaectively.  We  have  selected  u  =  0,  that  is,  two  dimensions,  and 
are  approximating  the  ideal  pattern  for  B  =  4  shown  in  figure  13.  It  is  seen 
that  the  approximations  become  progressively  better  as  e  decreases,  and  that 
the  result  in  figure  F-8  is  indeed  very  close  to  figure  13. 


The  approximation  w£(  s)  in  (F-12)  and  figure  F-5  used,  for  the  impulse 
area,  the  value  C  given  by  (F-7)  as  a  limit  of  (F-6)  for  snail  e.  A  better 
approach  is  not  to  use  the  limiting  value,  but  to  use  the  actual  value  of  the 
pertinent  function,  since  we  would  like  good  approximations  for  moderate 
values  of  e,  not  just  very  small  e.  This  procedure  is  considered  in  detail  in 
[12J,  with  the  result 
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(F-14) 


The  patterns  for  this  approximation  are  depicted  in  figures  F-9  through  F-ll. 
The  result  in  figure  F-ll  for  e  =  .1  is  now  better  than  the  result  for 
e  -  .001  in  figure  F-8;  and  all  we  have  done  is  to  modify  the  area  of  the 
impulse  at  s  =  1.  The  result  for  e  =  .15  in  figure  F-10  indicates  a  modest 
change  from  the  ideal  and  would  be  acceptable  in  some  cases.  The  program  for 
the  pattern  evaluation  is  listed  at  the  end  of  this  appendix. 


Another  possibility  is  to  relocate  the  impulse  in  (1-e,  1)  to  best 
approximate  the  ideal  pattern  in  figure  13.  More  generally,  a  shaped  narrow 
oulse,  which  i  s  concentrated  toward  the  boundary  at  s  *  1,  could  be  used; 
these  possibilities  are  discussed  further  in  [12]. 


Some  additional  results  involving  delta  functions  and  Vessel  transforms 
are  given  here  in  appendix  H. 


Figure  F-5.  Final  Approximation  to  Weighting  (50) 
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PROGRAM  FOR  PATTERN  EVALUATION 

10  Epi  =  .l  !  Pattern  for  Weighting  'P-iT,' 

10  Hu=0 

30  Nu=- 1 . 5 

40  Be =4 

50  DIM  G < 0 :  2 4 O '■> 

do  COM  U,  Be  ,  Mu,  No,  Mill 

70  M2 1  =Mu*2+ 1 

SO  T=2*Eps-Eps  -2 

90  <Nu*i  >*FNIriuxnuvNu+l ,  Bc*SQR\T  •  » 

100  fi=0 

llO  B= 1 -Epi 

120  FOR  I u  =  0  TO  240 

130  U= I u*P I  -'40 

1 4  O  S  ■ 1 F  M  3 1  H1  *  F  H  S 1  B  1 *  •  5 

150  N=2 

160  H  =  ( B - 0  >* . 5 

170  F  = B  -  0  >  '  3 

130  Vo=9E99 

190  T =0 

200  FOR  K  =  1  TO  N-l  STEP  2 
210  T  =  T  +  FNS <fl+H*K > 

220  NEXT  K 
230  S=S+T 
240  v='.:st-r::'*F 

250  IF  0BS  t  V-Vo >  <  =0B3 1 V  >  *  1 E-4  THEN  31 O 

260  Vo=V 

270  N=N*2 

280  H=H*.5 

290  F=F*. 5 

300  GOTO  190 

310  G<  Iu>*Fl*4FNJmjiaiU'.  Mu,  U)+V  !  Voltage  Re  =  port 
320  PRINT  Iu,G<.  Iu  > 

330  NEXT  1 14 

340  PLOTTER  IS  "98720" 

350  LIMIT  25,175,35,245 
360  OUTPUT  705; "VS5" 

370  SCALE  0,240,-70,0 

380  GRID  40,10 

390  PENUP 

400  FOR  I u=0  TO  240 

410  PLOT  I  u,  20*LGT  \  ABS <  G<  I u > ✓G <  0 >  >  ) 

420  NEXT  I <4 

430  PENUP 

440  END 

450  • 

460  DEF  FNS'S  • 

470  COM  U , Be , Mu , Nu, M2  1 

480  T=1-S*S 

490  T  1  =FN  Jnuxnu <  Mu,  Ut-S  > 

500  T£“FN Inuxnu*  Nu, Be  *SQR(T>  > 

510  RETURN  S' M21*T  Nu*Tl*T2 
520  FNEND 

530  i 


FNJnu'nu  and  FNInuxnu  are  listed  in  Append i  D 
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Appendix  G 

EVALUATION  OF  A  8ESSEL  INTEGRAL  VIA  RECURSION 
The  integral  of  interest  is 


v(u,3)  =  f  dsK(u,s/-^  IV(B^7), 


(G-l ) 


where  the  kernel  K  is  the  Bessel  function  as  given  in  (11).  We  have,  via 
[5;  9.6.28], 


dB 


g  (u.B) 

V 


I 


ds  K(u,  s)  B 


=  3  g^_1(u,3). 

(This  relation  actually  holds  for  any  kernel  K,  not  just  (11).) 
(G-l), 

9V( u,0)  =0  if  v  >  0, 

we  have  the  integral  recursion 


(G-2) 

Since,  from 


B 

gv(u,B)  =  £  dt  t  gv_1(u,t) 
0 

We  already  know  the  starting  case  of 

g0(u,B)  =  J  ds  s(5)  J  (us)  IQ 
see  (23).  Substitution  in  (G-3)  yields 


for  v  >  0. 


(G-3) 


(G-4) 
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Appendix  H 

TWO  BESSEL  INTEGRALS  THAT  YIELD  GENERALIZED  FUNCTIONS 

The  starting  point  is  the  Hankel  transform  pair  in  (16)  and  (16).  If  we 
let  w(  s)  =  J(s-a)  (where  a  >  o)  in  (15),  we  get 


g(u)  =  a(~j  J  (au)  for  u  >  0  . 

Substituting  (H-l)  in  (16)  then  yields  the  useful  relation 
06 

(duuJ(su)  J(au)  =  t£(  s-a) . 

J  u  u  a 

0 

On  the  other  hand,  if  we  let  the  weighting  be  a  doublet. 


vd  s)  =  i^(s-a),  then  g(u) 


-u 


(au). 


The  inverse  relation  (16)  yields 


QO 

f  du 

o 


J  (  su)  «J  ,  ( au) . 
u  w-l 


However,  since 

f(s)£'(s-a)  =  [f(a)  +  f‘(a)(s-a)  +  ...^  S' (s-a) 
=  f(a)ps-a)  -  f'(a)£(  s-a). 


then 

J  du  U2  su)  J  _j(au)  =  -  j$'(  s-a)  +  ^  S(  s-a). 

0 

Equations  (H-2)  and  ( H -6 )  are  the  desired  results. 


(H-l) 


(H-2) 


(H-3) 


( M-4 ) 


( H— 5 ) 


(H-6) 
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